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Abstract 

We study magnetic dipole (Ml) transitions between two quarkonia in the framework of non- 
relativistic effective field theories of QCD. Relativistic corrections of relative order v 2 are investi¬ 
gated in a systematic fashion. Non-perturbative corrections due to color-octet effects are considered 
for the first time and shown to vanish at relative order v 2 . Exact, all order expressions for the 
relevant 1/m and 1/m 2 magnetic operators are derived. The results allow us to scrutinize sev¬ 
eral potential model claims. In particular, we show that QCD excludes both contributions to the 
anomalous magnetic moment of the quarkonium induced by low-energy fluctuations and contri¬ 
butions to the magnetic dipole operators of the type induced by a scalar potential. Finally, we 
apply our results to the transitions J/i/> —> r/ c 7, Y(15) —> 7^7, Y(25) —> rjb(2S) 7, T(2 S) —> 7^7, 
77 6 (25) —> Y(15)7, hb(lP) —> Xbo,i{lP) 3 an d Xb 2 (lP) —> hb(lP)3 by assuming these quarkonium 
states in the weak-coupling regime. Our analysis shows that the «//?/> —► r] c 3 width is consistent 
with a weak-coupling treatment of the charmonium ground state, while such a treatment for the 
hindered transition T(25) —> r/b 7 appears difficult to accommodate within the CLEO III upper 
limit. 

PACS numbers: 12.38.-t, 12.39.Hg, 13.25.Gv 
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I. INTRODUCTION 


The non-relativistic nature appears to be an essential ingredient to understand the dy¬ 
namics of heavy quarkonia. It has been established soon after the discovery of the J/-0 in 
1974 by many subsequent phenomenological studies on numerous observables of the cc and 
bb bound states. Hence, heavy quarkonium is characterized by the interplay among the sev¬ 
eral supposedly well-separated scales typical of a non-relativistic system: the heavy quark 
mass m, the inverse of the typical size of the quarkonium 1/r ~ mv and the binding energy 
E ~ mv 2 , where v <C 1 is the velocity of the heavy quark inside the quarkonium. Nowadays 
the effective field theory (EFT) approach has become the paradigm to disentangle problems 
with a hierarchy of well-separated scales. Two effective field theories, non-relativistic QCD 


nn 


, have been developed in the last 


(NRQCD) [I!, || and potential NRQCD (pNRQCD) 
decade. Applications of these two EFTs have led to a plethora of new results for several 
observables in quarkonium physics (for a review see ^|). 

Among the observables that haven’t yet been considered in an EFT framework, are radia¬ 
tive transition widths. T hey have been studied so far almost entirely within phenomenologi- 
d,lmlallO, u, lid,lid,lid,lid,lid,lid, is, 19] (a sum rule analysis is provided in 


cal models 

For a recent review we refer to Eichten’s contribution in 
can be found in 22]. Mostly, the models are based on a non-relativistic reduction of some 


0 ). 


2l|. A textbook presentation 


relativistic interaction assumed on a phenomenological basis. Eventually, a potential model 
coupled to electromagnetism is recovered. In this work, we will describe radiative transi¬ 
tions in the language of EFTs. In particular, we will employ pNRQCD to study radiative 
transitions in a model independent fashion. 

Two dominant single-photon-transition processes, namely electric dipole (El) and mag¬ 
netic dipole (Ml) transitions, are of considerable interest. Since, for reasons that will become 
clear in the following, Ml transitions are theoretically much cleaner than El transitions, we 
will restrict ourselves to Ml transitions in this work. 1 

The kinematics of a transition H —* H'^ in the rest frame of H, where H and H' are 
two quarkonia, is described in Fig. |T| In the non-relativistic limit, the Ml transition width 


1 However, El transitions are the most copiously observed, because their rates are enhanced by 1/v 2 with 
respect to the Ml case. We will report about El transitions elsewhere. 
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FIG. 1: Kinematics of the radiative transition H —► H'"f in the rest frame of the initial-state 
quarkonium H . Mh and Mw are the masses of the initial and final quarkonium, and fc 7 = |k| = 
{Mfj — M^ i ,)/(2Mh) is the energy of the emitted photon. 

between two S-wave states is given by 
r _ 4 2 ^7 

n 3 Si^>n n So 7 o o 

3 v rn z 

where eeg is the electrical charge of the heavy qnark (e& = —1/3, e c = 2/3), a is the 
fine structure constant and R n i(r ) are the radial Schrodinger wave functions. The photon 
energy /c 7 is about the difference between the masses of the two quarkonia, therefore, it is 
of order mv 2 or smaller. 2 Since r ~ 1 /(mu), we may expand the spherical Bessel function 
jo(k 1 r/2) = 1 — (/c 7 r) 2 /24 +.... At leading order in the multipole expansion, for n = n', the 
overlap integral is 1. Such transitions are usually referred to as allowed. At leading order, for 
n ^ n the overlap integral is 0. These transitions are usually referred to as hindered. The 
widths of hindered transitions are entirely given by higher-order and relativistic corrections. 

Equation 0 is not sufficient to explain the observed transition widths. In the case of 
allowed ones, for instance, it overpredicts the observed J/^ —> ?/ c 7 transition rate by a 
factor 2 to 3. A large anomalous magnetic moment or large relativistic corrections have 
been advocated as a solution to this problem. Hence, it is crucial to supplement Eq. 0 
with higher-order corrections. EFTs provide a systematic and controlled way for doing it. 

2 This is in sharp contrast with radiative transitions from a heavy quarkonium to a light meson, such as 
J/V’ —»777, whereas a hard photon is emitted. 


drr 2 R n , 0 (r) R n 0 (r) j 0 


k 7 r 


( 1 ) 


3 








EFTs are characterized by a power counting and a range of validity (the system must 
consist of a specific hierarchy of scales). Errors are controlled by the power counting; higher- 
order corrections can be systematically included. Among these, EFTs include corrections 
coming from higher-Fock states, typically missed in potential models. In particular, in both 
NRQCD and pNRQCD color-octet contributions play a crucial role in some processes. 

NRQCD is obtained from QCD by integrating out modes of energy m. The energy 
scale m is sometimes called hard. We denote with Aqcd the typical hadronic scale. Since 
m 3 > Aq CD , the matching procedure that ensures the equivalence of the two theories may 
be carried out in perturbation theory. At this stage, also hard photons are integrated out. 
However, at the accuracy we are interested in, their contribution is negligible. 

pNRQCD is obtained from NRQCD by integrating out modes of energy mv. This scale 
is sometimes called soft. We shall distinguish between strongly coupled quarkonia, for which 
mv ~ Aq CD and weakly coupled quarkonia, for which mv 2 > Aqcd- bi the first case, the 
matching has to be done in a non-pertnrbative fashion. In the second case, it may be done 
order by order in the strong-coupling constant. Low-lying quarkonia are believed to be in 
the weak-coupling regime, higher excitations in the strong-coupling one. Soft photons are 
also integrated out at this stage, but its contribution is numerically irrelevant with respect 
to that one coming from soft gluons. In the strong-coupling regime, the degrees of freedom 
of pNRQCD (coupled to electromagnetism) are singlet quarkonium fields and photons of 
energy and momentum of order mv 2 or smaller. The scale mv 2 is sometimes called ultrasoft. 
In the weak-coupling regime, there are also octet quarkonium fields and ultrasoft gluons. 
Ultrasoft fields are multipole expanded about the centre-of-mass coordinate. The power 
counting of the pNRQCD Lagrangian goes as follows. Ultrasoft gluons and virtual photons 
scale like mv 2 , the real photon, emitted in a single-photon transition, scales like mv 2 or 
smaller. In addition, the matching coefficients inherited from NRQCD are series in a s . To 
simplify the counting, we will assume that a s {m) ~ v 2 . In the weak-coupling regime, the 
matching coefficients of pNRQCD can be calculated in perturbation theory. Since the static 
potential is proportional to a s (l/r)/r ~ mv 2 , it follows that a s (l/r) ~ v. 

In the paper, we will mainly work out pNRQCD in the weak-coupling regime. Therefore, 
our final expressions will be applicable only to the lowest quarkonium resonances. However, 
some intermediate results will also apply to the strong-coupling regime. In particular, the 
1/m and 1/m 2 matching will be valid to all orders in a s . 
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Some of the results presented here are new, some may be understood as a rewriting in the 
language of EFTs of results already derived long time ago in the framework of phenomeno¬ 
logical models. Among others, we will address and answer the following questions. (%) What 
is the size of the quarkonium anomalous magnetic moment? (ii) Is there a scalar interac¬ 
tion contribution to Ml transitions? (in) What is the size of the octet contributions to Ml 
transitions? We will end up with a rather concise formula which takes into account the full 
0(k 3 v 2 /m 2 ) relativistic corrections. We will clarify the validity and range of applicability of 
the widely-used formula of Ref. ^|. Applications to some Ml transitions between low-lying 
quarkonia will be discussed at the end. 

The paper is organized as follows. In Sec. El we first briefly review NRQCD and pN- 
RQCD, then work out the basic formalism and calculate the transition widths in the non- 
relativistic limit. In Sec. m we match the electromagnetic interaction Lagrangian of pN- 
RQCD relevant for Ml transitions up to 1/m 3 terms. In Sec. IIVI we calculate contributions 
to the transition widths from wave-function corrections and, in particular, color-octet con¬ 
tributions. In Sec. El we sum all corrections and give the final formulae valid up to order 
fc 3 v 2 /m 2 . In Sec. ED the decay rates of J/-0 —■» ij c 7 , T(IS') —■» r] b 7 , T(2S') —> r) b (2S)'y, 
T(2S') -> 77 b 7 , rjh(2S) -> TITS') 7 , h b (lP) -> x 60 i i( 1 P )7 and y fe2 ( 1-P) -»• h(lP)^ are cal¬ 
culated. Finally, in Sec IVIII we conclude. In one appendix, we discuss alternative ways to 
derive final-state recoil effects, in the other one, issues about gauge invariance. 


II. MAGNETIC DIPOLE TRANSITIONS: BASIC FORMALISM 

A. NRQCD 

NRQCD is the EFT that follows from QCD by integrating out hard modes, i.e. modes of 
energy or momentum of order m fl, 2], To describe electromagnetic transitions, we need to 
couple NRQCD to electromagnetism. For simplicity, we will call this new EFT also NRQCD. 
The effective Lagrangian is made of operators invariant under the SU(3) C x f/(l) em gauge 
group. We display here only the part of the Lagrangian, which is relevant to describe Ml 
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transitions at order k^v 2 /m 2 \ 


C 


NRQCD 


= ^ ( iD ° + L) v ’ + 


a ' 9B ^ ~ 8m 2 ^ (T ' H Dx > 0 E ]^ 
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■ ee Q B em 0 - ~^i' ] cr ■ [-iDx, ee Q E em ]0 
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8 m 3 


ibU D 2 , rr ■ ee o B em }0 - ^DV • eeoB cm D> 

4 m 3 


JOL,l, t 


0 T (<r • D ee 0 B cm D + D ee Q B em a ■ D) 0 


8 m 3 
+ [0 -> 


iT /tem 
*/i 5 




“1“flight ; 


( 2 ) 


where 

Aw = - jF' “‘F^, - + £ 5/ If . (3) 

/ 

and 0 is the Pauli spinor held that annihilates a heavy qnark of mass m, havor Q and 
electrical charge eeg, x is tl ie corresponding one that creates a heavy antiquark, and q/ are 
the light quark Dirac fields. The gauge fields with superscript “em” are the electromagnetic 
holds, the others are gluon helds, iD 0 = id 0 — gT a A g — eeQA™, ?D = i V + gT a A a + ee Q A em , 
[Dx.E] = D x E - E x D, E' = F i0 , B' = -e ijk Fi k / 2 , E iem = F i0em and B iem = 
-<i jk FA'""/2 (e 123 = 1 ). 

The coefficients cp, eg, cfd, cf 11 , c^\, cf^ 2 and c®™ are the matching coefficients of the 
EFT. They satisfy some exact relations dictated by reparameterization (or Poincare) invari¬ 
ance 1231: 


c e s m = 2 c e F m - 1 , c s = 2 c F - 1 , 

,,em _ em _ i 

L W 2 — C W1 1 5 

r em _ em i 

c p'p ~ C F 1 - 


(4) 

(5) 

( 6 ) 


Note that the c| 0 ?: are independent of cfP. All the coefficients are known at least at one loop 
2^ |. In particular, we have 

cT = 1 + 4 m = 1 + + OK 2 ), (7) 

c - = 1 + c ^(^ + l h v) + 0( “- 2) ’ (8) 

where Cp = (iV 2 — l)/(2A r c ) = 4/3 and ( 7,4 = N c — 3; /tg 11 is usually identified with the 
anomalous magnetic moment of the heavy quark. Since and c“ n are 1 + C?(a s ), c ^ 2 
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and Cy£ are 0(a s ). Kq 11 is less than 10% for charm and bottom. One may expect that the 
magnetic moment of the quarkoninm may be larger than that, because, apart from inheriting 
the magnetic moments of the quarks, it may get potentially large low-energy contributions. 
We will clarify this point in Sec. 1111 Bl 

In general, the matching coefficients will contain contributions coming from virtual pho¬ 
tons of energy or momentum of order m, which have also been integrated out. These 
contributions are suppressed by powers of a and shall be neglected in the following. We will 
only consider QCD corrections. 

B. pNRQCD 

NRQCD still contains redundant degrees of freedom in describing a quarkoninm state 
far below the open flavor threshold. pNRQCD is the EFT that follows from NRQCD by 
further integrating out quarks and gluons of momentum and energy of order mv and gluons 
of momentum of order mv and energy of order mv 2 jjj 0]. We consider here pNRQCD in 
the weak-coupling regime (mv 2 > Aq CD ). The degrees of freedom of pNRQCD are quarks 
of momentum mv and energy mv 2 and (ultrasoft) gluons of energy and momentum of order 
mv 2 . Since we are interested in quarkonium states, it is convenient to express the pNRQCD 
Lagrangian in terms of quark-antiquark fields. These are 3 ® 3 tensors in color space and 
2 ® 2 tensors in spin space, which depend on the centre-of-mass coordinate R and the 
relative distance r of the two quarks. At leading order, the pNRQCD Lagrangian very much 
resembles a potential model, where the potentials are the matching coefficients of the EFT 
that encode the soft-scale contributions. The pNRQCD Lagrangian, however, also contains 
dynamical (ultrasoft) gluons and their interactions with the quark-antiquark fields. 

Quarkonium radiative transitions involve real photons of energy and momentum of 
order mv 2 for hindered transitions and smaller for allowed ones. These transitions are 
described by pNRQCD if photons of momentum mv are integrated out from NRQCD and 
photons of energy and momentum of order mv 2 or lower are explicitly coupled to the quark 
fields in the pNRQCD Lagrangian. To ensure that gluons and photons are of energy and 
momentum not larger than mv 2 , all gauge fields are multipole expanded in the relative 
distance r, and, therefore, depend on the centre-of-mass coordinate R only. 

Gauge invariance can be made manifest at the Lagrangian level by reexpressing the quark- 
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antiquark fields in terms of fields that transform like singlets under f/(l) em and like singlets 
or octets under SU (3) c gauge transformations. We denote these fields as S = S \J a J~N~ c and 
O = \[20 a T a , respectively. 

The pNRQCD Lagrangian, which is relevant to describe Ml transitions at order k 3 v 2 /m 2 , 
is given by 


C 


pNRQCD 


d 3 r Tr <j S f (id Q + ^ - V s ) S 

[ \ 4m m ) 

f D 2 V 2 \ 

+0 t (w 0 + — + - v 0 o 

\ 4 m m ) 

+V A (O f r ■ gE S + S f r • gE O) + V B ( Qt ’ r '^ E } q 


+£ 7 pNRQCD 
+Aight ; 


(9) 


where 


C 


7pNRQCD 


d 3 r Tr 1 VT S f r • ee Q E cm S 


{S\*-ee Q B em }S 


1 cr-B 

- 

2 m 5 

1 tJ et * j 

ItaW b''”' 


r'r'i V'V'rr • ee Q E cm )} S 
{0*. CT . eeo B'”}0 

<r-(r Xr X B) 

1 Vo m a 


4 m 2 


1 VJ 


CTvB 


{SU'[rx(rx eegB em )]} S 
{S^, <r • eegB om } S 


4 m 2 r 

1 ( 7 -VxE 

h^'HVx.eeoEyJS 

1 cr-'VrX'r-'V E 

Vc ^ 


16 m 2 5 

1 zi 

J_T/ ™ 

4 m 3 51 

| (Vr-cr) (Vr-B) 

j _T/ m 3 

+ 4m 3 ^ 


[S*,[-iV^r^eegE^JS 
{St, <r • eegB 6 ™} V 2 S 
{St, cr* eenB emj ] 


If not differently specihed, all gauge fields are calculated in the centre-of-mass coordinate 
R, iD 0 O = id 0 O - g[T a A%, O], iDO = iVO + g[T a A°,0], V* = d/dK 1 and V* = d/dv 1 . 
The trace is over color and spin indices. 



In the initial quarkonium rest frame, the power counting goes as follows: 3 V r ~ mv , 
r ~ 1 /mv and E. B ~ m 2 v A . The electromagnetic fields associated to the external photon 
scale like E cm , B cm ~ k 2 . The centre-of-mass derivative V acting on the recoiling final 
quarkonium state or emitted photon is of order k T Operators that have not been displayed 
are suppressed either in the power counting (e.g. 1/m 4 singlet operators) or in the matching 
coefficients (e.g. a l/(m 3 r 2 ){S^ a • eegB cm }S operator does not show up at tree level) or 
because they project on higher-order Fock states (e.g. 1/m 2 octet operators). 

The coefficients V in Eqs. © and m are the matching coefficients of pNRQCD. The 
matching coefficients of Eq. © have been calculated in the past years. We refer the reader 
to P| and references therein. In the following, we will calculate the matching coefficients of 
Eq. (g5D . Here, we only note that since mv Aq CD , they may be calculated in perturbation 
theory. Vs and Vo play the role of a singlet and octet potential. They may be arranged in 
powers of 1/m. The static contribution is the Coulomb potential: 


Vf = ~C F ^. 


E (0) - 

“2 Nr r 1 


1 a Vo 


( 11 ) 


where, at leading order, ay s = ay 0 = a s . In a Coulombic system a s (l /r) ~ v. 

Let us discuss the different terms appearing in Eqs. © and (hod- The first three lines 
of Eq. © display the pNRQCD Lagrangian in the limit of zero coupling to the photons. 
The third line describes the coupling of the quarkonium fields to ultrasoft gluons at order 
r in the multipole expansion. Higher-order terms are irrelevant for the present purposes. 
At tree level, the coefficients Va and Vb are equal to 1. Equation m provides the part of 
the pNRQCD interaction Lagrangian coupled with the electromagnetic held relevant for Ml 
transitions. The first line describes the coupling of the quarkonium singlet held to ultrasoft 
photons at order r. This is the familiar El transition operator (Rf™ = 1 at tree level). As 
we will discuss in the following, if the recoiling of the hnal-state quarkonium is taken into 
account, this term contributes to Ml transitions. A similar term involving the coupling 
with the octet held is suppressed in the transition amplitude. From the second line on, we 
display spin-dependent operators coupled to ultrasoft photons. Lines two and three come 
from multipole expanding the magnetic dipole operator at O(r 0 ) and 0(r 2 ) respectively. 
The order r term does not contribute to Ml transitions. Line four represents the leading 
magnetic dipole operator for octet quarkonium helds. 


3 For simplicity, we give the power counting in the case mv 2 ~ Aqcd only. 
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C. Radiative transitions 


The process H —a H' 7 is described in the rest frame of the initial quarkonium state H 
by the kinematics of Fig. |T| The transition width is given by: 


H^H' 7 


J3p/ 1 

( 2^)3 ( 2^)3 2 k^ S ^ PH ~ k ~ P ) 


X 




El-4 [#(<>, A) ^ff'(P',A'b(k,<T)] 


Wo 


1 


87T 2 


Mj. 


H J J 0 


dkk dQ(k)8(k-k^ 


x 


^ E M (0. A ) - A'h(k, 0-)] I 2 . (12) 


Wo 


where = (^/P /2 + Mjj,, P'), AT = (|k|,k) and 
-4 [P(0 , A) -a P'(-k, A')7(k, a)] (27t) 3 5 3 (P / + k) 

= —(H'(P', A 7 )7(k, cr)| J d 3 P£ 7pNRQCD |^(0 ,A)). (13) 

In Eq. m , the initial state is averaged over the polarizations, whose number is N\. 

The quarkonium state \H(P. A)) is an eigenstate of the pNRQCD Hamiltonian with the 
quantum numbers of the quarkonium H. It has the non-relativistic normalization: 


(H( P', A')|P(P. A)) = Sw (27r) 3 <5 3 (P - P'). 

The photon state I'y(k, cr)) is normalized in the usual Lorentz-invariant way: 

(7(k, cr)|7(k / ,cr')) = 2k 8 m' (27r) 3 5 3 (k - k'). 


(14) 


(15) 


D. Quarkonium states 


According to the power counting, the leading-order pNRQCD Hamiltonian is given by 


H. 


(o) 


pNRQCD 


= [ d?R j d 3 r Tr {s* hf S + 0 f fcg* o} + H llgU , 


(16) 


where 


hf = -^+y o) 

m 


s 1 


*,<») = -^ + y°> 

m 


o 1 


(17) 
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and //light is the Hamiltonian that corresponds to £u g ht- The spectrum of pNRQCD has 
been first studied in jj], to which we refer for discussions. We call |//(P, A))® the subset of 
eigenstates made by a quark-antiquark pair in a singlet representation: 


|//(P, A)) (0) = J d 3 R J d 3 r e iPR Tr|<f> 


(o) , 
HW 


S f (r,R)|0)} , 


(18) 


where |0) is a state that belongs to the Fock subspace containing no heavy quarks but an 
arbitrary number of ultrasoft gluons, photons and light quarks. The state |0) is normalized 


in such a way that Eq. m is fulfilled. The function $ 


(o) , 

H{ A) 1 


= <0|S(r,R)|ff(0,A)><°> is an 


eigenstate of the spin and orbital angular momentum of the quarkonium and satisfies the 
Schrodinger equation 

hf = 4 0) *2kW. (19) 

is the leading-order binding energy of the quarkonium H\ Mh = 2 m + E^. For later 


use, we write for L = 0 states, 


do) 




1 R n o(r)' T ' e ”’ Sl(A) 


s/2 


♦SkW = 


1 


7S RM{r) s/r 


( 20 ) 

( 21 ) 


where e n 3 Sl (A) is the polarization vector of the state n 3 Si, normalized as e* 3Si (A) -e n 3 Sl (A 7 ) = 
S\y, and for L = 1 states, 


$ nip l( A)( r ) = 


P 0 ( r ) = 



3_ e n i Pl (\)-r 

in Knl[r) ^2 

, 1 a ■ r 


d> 


XW r ) = \!^ R ^ 


<J)( 0 ) (V) — 

n 3 P 2 ( A)f > 


i/ Rnl{r) 


<r • (r x e n 3p 1 (A)) 

y/2 

(A)fl 

s/2 ' 


( 22 ) 

(23) 

(24) 

(25) 


where e n i Pl (A) and e n 3 Pl (A) are polarization vectors satisfying e* iPi (A) -e„i 

^72 3 Pl (A)- 

e n 3 Pl (A 7 ) = 5\\i, whereas the polarization of the n 3 P 2 state is represented by a symmetric 
and traceless rank-2 tensor hR 3p2 ( A) normalized according to /q) 3 *p 2 (A) h 3 n 3p2 (A 7 ) = 5\\>. 

The state |//(P, A)) may be obtained from |//(P. A))^ by quantum-mechanical pertur¬ 
bation theory. At relative order v 2 the following corrections may be relevant. 
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1. Higher-order potentials 


In the weak-coupling regime, corrections to the zeroth-order pNRQCD Hamiltonian of 
the type 

5H = j d 3 R j d 3 r Tr {SUl^S} , 

are typically suppressed by v 2 with respect to the leading term. First-order corrections to 
the quarkonium state, induced by these terms, are therefore relevant at relative order v 2 : 

|#(P,A)>« = [ d 3 R [ d 3 r e* PR Tr{^ (A) (r)St(r,R)|0)} , (26) 


where 

Cj-) 

s*h(>.){t)= Y. rwrfim (h\\')\svs\h(\)), ( 27 ) 

E H i 

and (r\H(X)) = $® A )(r). Here and in the following, we shall use the Dirac ket to indicate 
the eigenstate either of a quantum-mechanical operator (like |r), which stands for an eigen¬ 
state of the position operator, or \H(X)), sometimes also written as \nL), which stands for 
an eigenstate of h^) or of a quantum-field operator (like |ff(P,A)), which stands for an 
eigenstate of the pNRQCD Hamiltonian). 

In general, 8V$ may also depend on the centre-of-mass momentum P. We shall distin¬ 
guish between zero-recoil corrections (where SVs does not depend on P) and (final-state) 
recoil corrections (where SVs depends 011 P). Effects of these corrections to the transition 
amplitude will be discussed in Sec. IfV Al 


2. Higher Fock-space components 

The leading correction to the quarkonium state that accounts for the octet component is 
induced by 

SH=- J d 3 R J d 3 r Tr {O f r • gE S + S f r • gE 0} . (28) 

According to the power counting, this is a correction of relative order v. The first-order 
correction to the quarkonium state is 

|AT(P, A)> (1) = f d 3 R f d 3 r e iFR ! d 3 x Tr J O f (r, R)(r| —-^-|x) 

J J J ^ Hq -Hlight 

x[-x-sE(R)] *g> (J) (x)|0)|. (29) 
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Since it has a vanishing projection on |i7(P, A))(°q in a transition matrix element it con¬ 
tributes at relative order v 2 . Second-order corrections are of relative order v 2 . They contain 
two orthogonal parts: 


I#(P, A)) <2 > = |77(P, A))f + \H(P, A))* 2 *, (30) 


where 


|tf(P.A))f = d 3 R d 3 r e iP R 


d 3 x 
(o) mW°) 


J d 3 y Tr -f S f (r, R 


x ( r i E JJ y 1: E V > [-y&m 


x(y| 


H'^H, A' 

1 


771(0) pi(0) 77 

-C/77 — &H' — -flight 


Ix) [—x • f?E(R)] $2(A)(x) |0) , 


^(°) _ 7,(0) 

xj h n o 


H, 


light 


(31) 


\H(P, A)){, 2) = \H(P, A))(°), (32) 

^1 + E|E) ( 2 vr )353 ( p_p/ )= (°)(#( P ' ; A)|ff(P, A)) = (2 tt) 3 5 3 (P - P') .(33) 

l+hZ//(A) is the usual quantum-mechanical normalization constant of the state, which stands 
for the probability to find the leading color-singlet component in a physical quarkonium state. 
Effects of these corrections to the transition amplitude will be discussed in Sec. IIVBI 


E. Ml transitions in the non-relativistic limit 

In accordance with the power counting of pNRQCD, the leading contribution to Ml 
transitions comes from 

4!?i = Jd 3 r TV {2- V'f {St, <r - eeo B“} S j . (34) 

a-B 

As we will discuss in the next section, at leading order V s m =1. 

For 5-wave quarkonium, substituting Eqs. CEBD, EDD, (EH) and (EH into Eq. m leads 
to 

M (0) [n 3 5i(0, A) -v n' 1 S’ 0 (-k) 7 (k, cr)] = 8 nn > iee Q Gn3Sl ^ (k x e (a)) , ( 35 ) 

L J m 
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where we have used that 


( 7 (k,a)|B em (R)|0) = -ikx e*(cr) e _ik R . (36) 

poo 

The factor 5 nn > in Eq. (1351) comes from the overlap integral / dr r 2 R n e(r) R n ^(r). Substi- 

J o 

tuting the transition amplitude into Eq. m, we obtain: 


rn 3 Si—>n' 1 S'o 7 8 nn ' Oi6, 


US 

2 K "f 

3 Q m 2 


1 - 




M, 


n 3 S i 


(37) 


The term — 5 nn / k 1 /M n 3s 1 is negligible at order k^v 2 /m 2 : it vanishes for hindered transitions 
and is of order v 4 for allowed ones. Hence, Eq. (1371) gives back Eq. m at leading order in 
the multipole expansion, i.e. the well-known formula of the transition width in the non- 
relativistic limit. 

For P-wave quarkonium we obtain: 

e*,i Pl (A') • (k x e*(cr)) 


«4 (0) [n 3 P 0 (0) —> n n P 1 (-k, A') 7 (k, a)\ = 8 nn 'iee Q 
*4 (0) [n 3 Pi(0, A) - n n P 1 (-k, A') 7 (k, cr)] 

8nn' 




m 


e*/i Pl (A') • [e n 3 Pl (A) x (k x e*(<r))] 


Al (0) [n 3 P 2 (0, A) - u ,1 P 1 (—k, A') 7 (k, a)] 

A • e; 
Onn' 




m 


)n Pl (A') hP 3p2 (X) (k x e*(cr)) J 


m 


(38) 


(39) 


(40) 


Substituting the transition amplitudes into Eq. ED, we end up with 


^-'n 3 Pj—>n’ x Pi 7 8 nn i 01 €q ^ k^ ( 1 


Ag 


(41) 


M n 3 Pj _ 

The width for P-wave spin-singlet to spin-triplet transitions is obtained by multiplying the 
right-hand side of Eq. (l-TT!) by (2 J + l)/3. 

In the following, we will concentrate on higher-order corrections to 5-wave transitions. 
We shall come back to P-wave transitions in Sec. IYD 


III. MATCHING OF PNRQCD MAGNETIC DIPOLE OPERATORS 

Our aim is to complete Eq. ED with corrections of relative order v 2 . In accordance with 
the counting a s {m) ~ v 2 and cx s (l/r) ~ v, these also include corrections to the matching 
coefficients of pNRQCD. 
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FIG. 2: Green function with four external quark/antiquark lines and a photon. The time flows 
from the bottom to the top. 


The matching coefficients of pNRQCD encode glnons of energy or momentum of order 
mv. Since this scale is associated with the distance between the two heavy quarks, the 
matching coefficients are, in general, functions of r. They also contain hard contributions, 
typically encoded in the matching coefficients inherited from NRQCD. In the following, we 
will retain the full matching coefficients of NRQCD and count the matching from NRQCD 
to pNRQCD only in powers of a s calculated at the soft scale. We will exploit the explicit 
form of the NRQCD matching coefficients in Secs. 0 and I VI I 

The matching from NRQCD to pNRQCD may be performed by calculating Green func¬ 
tions in the two theories and imposing that they are equal order by order in the inverse of the 
mass and in the multipole expansion. In particular, to match the electromagnetic couplings, 
we need Green functions with four external quark/antiquark lines and one external photon 
line as shown in Fig. El The matching condition reads: 



(42) 


Since we are working in a situation where the typical momentum transfer between the heavy 
quarks is larger than Aqcd ; we may, in addition, perform the matching order by order in a s . 

A. Matching at 0(1) 

If we aim at calculating the matching at 0(1) a convenient approach consists in projecting 
NRQCD on the two-quark Fock space spanned by 



( 43 ) 
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where ip{x\,x 2 ,t) is a 3 0 3 tensor in color space and a 2 0 2 tensor in spin space. After 
projection, all gluon fields are multipole expanded in r. Gauge invariance is made explicit 
at the Lagrangian level by decomposing 


<p(xi,x 2 ,t) 

S'(R,r,t) 

0’(R,r,t) 


P exp ^ig j A ■ dx^j S'(R, r, t) 

+P exp ^ ig j A • dx^j O^R, r, t) exp ^ ' ig J A ■ dx^j , 
exp ^ieeQ J A em • dx^j 5(R, r, t), 

exp (ieeQ f A em • dx] 0(R, r,t), 


(44) 

(45) 

(46) 


\ Jxi / 

where P stands for path ordering, R = (xi + X2)/2 and r = xi — X 2 . The helds S and 
O transform like singlets under f7(l) em gauge transformations and like singlet and octet 
respectively under SU(3) C gauge transformations. After projecting (J2J) on (TIH1) . one obtains 


vr = i, 

(47) 

a-B 

T/ m _ /m 

V S — C F ? 

(48) 

( r cr-B 

w ’ rn _ „em 

V S — c F 1 

(49) 

cr-B 

T/ m _ ,em 
y O ~ C F ? 

(50) 

cr-VxE 

T/ m 2 _ „em 

V S — C S > 

(51) 

cr-Vr-Xr-VE 

T/ ^ 

V S — C S > 

(52) 

a-B 

T/ m3 _ em em _ i 

V S — C W1 L W 2 — 1 ) 

(53) 

(V r -<t) (V r -R) 
t/ ^ _ ,,em 

— Vp ' 

(54) 


cr-(rXrxB) cr-B 

The matching coefficients V s m and V s m are zero at 0(1). 

We consider now the impact of the 0(1) matching on the transition amplitude. In order 
to keep the notation compact, it is useful to define 


A = 


A 

A® ’ 


(55) 


where A is an amplitude calculated from Eq. (THT1) by substituting £ 7P nrqcd with the con¬ 
sidered operator and A*' 0 ' is connected with the leading Ml amplitude M (0) (|SS1) by 


S„„,A ,0) = -4< 0 >. 


(56) 
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(A) The matching coefficient (1481) induces the following correction to 


A [VS^O, A) -> n , 1 S' 0 (-k) 7 (k,cr)] = k™ 5 n 


( 57 ) 


(B) The correction induced by the operator ° F {S^, r'r^ V l V J cr • eenB em )| S 

16m 


is 


A [n 3 S'i(0, A) —> n / 1 S' 0 (-k) 7 (k, a)] = — k 2 (n'S\r 2 \nS) 


(58) 


c 

(C) The correction induced by the operator — 5 2 [S^, cr • [— i'V x, ee^E 6 ™]] S is 


A [n 3 5i(0, A) -> n' 1 S' 0 (—k)y(k, cr)] = c e s m — 8 n ' n , 
L J 8 m 


(59) 


where we have used 

(y(k, a)|E cm (R)|0) = -ike*(a)e~ ik R . (60) 

c em 

(D) The correction induced by the operator-TS 1 ”, cr • T—fV r x, r*(V*eenE em )ll S is 

16 m z 

A [n^O, A) —> ?/ 1 S' 0 (-k) 7 (k, a)] = 4“ ^ ^ 5 n f n + (rAS^r • plrnS 1 ) ^ • (61) 

(E) The correction induced by the operator-- {S^. cr ■ eenB cm } VjiS is 

Am 6 

A [n 3 5i(0, A) — ► n / 1 5 ' 0 (-k) 7 (k,(j)] = ~{n'S\-^-^\nS) . (62) 

c em 

(F) The correction induced by the operator {st, cR eegB 6 ” 1 - 7 } V* V^S is 

A [n 3 5i(0, A) -*• n , 1 S' 0 (-k) 7 (k,a)] = — ^-(n'S\^-^\nS), (63) 

where we have used c= k'q 1 . 


cr-B 

B. Calculation of V s m 

In this section, we match the operator (EH beyond 0(1). This operator provides the lead¬ 
ing transition widths dSZD and m in the case of allowed Ml transitions. Hence, corrections 

a-B 

of order a s and a 2 to V s rn , which may arise in the matching from NRQCD to pNRQCD, 
are potentially larger than or of the same order as genuine relativistic v 2 corrections to the 
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transition width. Surprisingly, we shall be able to perform the matching exactly and provide 
a result that is valid to all orders in perturbation theory and non-perturbatively. 

Before going to the matching, we recall that the matching coefficient c e ™ that appears in 
the NRQCD Lagrangian © is the heavy-quark magnetic moment. The matching coefficient 

a-B 

V s m is the magnetic moment of the singlet-quarkonium held. While the first one gets only 
contributions from the hard modes, the second one may potentially get large contributions 
steming from the soft scale. 

a-B 

The matching of V s m proceeds as follows. 

(1) First, we note that the only amplitudes in NRQCD that may contribute to the 
matching are those where the photon couples to the heavy quark (antiquark) through the 
operators tfya ■ eegB cm ^/m or — x*a ■ eeQB cm x/m. At order 1/m, this is the only magnetic 
spin-hipping coupling to the quark. If the photon couples to loops of massless quarks then, 
at leading order in the electromagnetic coupling constant, the sum of the electric charges 
over three light havors vanishes. In the bottomonium system, contributions from charm- 
quark loops should also be considered. If the momentum flowing through the loop is hard, 
the contribution is suppressed by, at least, a^mb) and, therefore, beyond the accuracy of 


this work. If the momentum is soft, the charm quark effectively decoup 
can be described by an effective held theory with three massless quarks 


es and the system 


MM 


(2) The crucial point is to recognize that there are no extra momentum or spin dependent 

a-B 

insertions on the heavy-quark lines that contribute to the matching of V s m , since they carry 
extra 1/m suppressions. As a consequence, the magnetic spin-hipping operator cr-eeQB om /m 
behaves, for the purpose of the matching, as the identity operator in coordinate space 4 and 
the magnetic matrix element factorizes. Therefore, any (normalized) NRQCD amplitude 

a-B 

contributing to V s m may be written as 

em rtf 


2 m 


dt ( 7 |cr U) ■ eeQB em (x!,£) + fr {2) ■ eeQB em (x 2 ,£)|0), 


(64) 


which, after multipole expansion of the magnetic held, becomes (we neglect terms propor- 


4 The property that <x • eeQB em /m depends neither on gluon fields nor on the relative coordinate r, which 
in turn is a consequence of the ultrasoft nature of the external photon, will be used again and again in 
the course of the paper and is responsible for most of results. 
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r tf 


dt 'X/X/X/X/t] t 


em cee Q B c 
C F 


1 1 


t2 


2 m 


'X/X/X/X/'H 


+ 


+ 


'X/X/X/X/t] 


X 



'x/x/x/x/TD 


FIG. 3: Diagrammatic factorization of the magnetic dipole coupling in static NRQCD: tf > t\ > 
t 2 > ti. Dashed lines are longitudinal gluons, the box represents the magnetic dipole coupling. 

tional to r): 

( em em \ rtf 

h + jL (r ' V)2 + ''') (' r “ ) + °’ (2> ) ' l *<7ke Q B'“(i)|0>, (65) 

where <r (1 ) stands for the Pauli matrices acting on the quark and (t (2> for the Pauli matrices 
acting on the antiquark. 

To see how factorization works from a diagrammatic point of view, let us consider a 
photon insertion on a quark line. In general, this happens in between two longitudinal gluon 
insertions (first diagram of Fig. |2J). Transverse gluons couple to quark lines through 1/m 
suppressed operators and are irrelevant for the purpose of the matching. At order l/m°, the 
coupling of longitudinal gluons to quark lines is spin and momentum independent. Therefore, 
cr ■ ee^B 6111 may be freely moved along the quark lines. The sum of the first three diagrams 
of Fig. 01 is proportional to 

0(tf - ti)9(ti - t)9(t - t 2 )9{t 2 - ti) + 9{tf - t)9(t - t 1 )9(ti - t 2 )9{t 2 - U ) 

+9{t f - t 1 )9(t 1 - t 2 )9(t 2 - t)9{t - ti) 

= 9{t f - t)9(t - ti)9(t f - ti)9(ti - t 2 )9(t 2 - ti ), 
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(a) 




fb) 


D 2 

2m 


\ 


cf- 


cr_gB 

2m 



cr-(rXrxB) 

FIG. 4: Diagrams contributing at order a s to F s “ . Although not displayed, the symmetric 

diagrams are understood. D is the covariant derivative under SU( 3) c x C/(l) em - 

where the thetas come from the static heavy-quark propagators. The equality completes the 
factorization proof graphically represented in Fig. |oJ 
(3) By matching Eq. (1(751) to the pNRQCD amplitude 


2 m 

it follows that 


£_£ ( r .y)2om 

V s m m , 2 

5,5 f v . _|_ 


+ 


16m 


(< 


r(l) 


r(2) 


)•/ dt (7|ee Q B em (f)|0), (66) 


V s m =Vg ’ m =4 m , 


(67) 


Equation (HD is a result that holds to all orders in the strong-coupling constant and also 
non-perturbatively. It excludes that the 1/m magnetic coupling of the quarkonium field is 
affected by any soft contribution. A fortiori, it excludes large anomalous non-perturbative 
corrections to this coupling. 


cr-(rXrxB) a_ 

Calculation of V s m and Vg 


q(rXrxB) 

The matching coefficients V s m and Vg 


' S ami v s do not get contributions at 0(1). At order 

a s , the diagrams contributing to the NRQCD part of the matching are shown in FigQ] Note 
that the photon is emitted by the electromagnetic field embedded in the covariant derivative. 
If we sandwich the diagrams between initial and final states that are gauge invariant under 
SU( 3) c x I/(l) em and multipole expand the external electromagnetic field, we obtain: 

cr(rXrXB) 


Vo 


CfOL S , , CfOIs 

—— (2 c F - c s ) = —— 


V s m = 0. 


2 
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( 68 ) 

(69) 


















In the first equation, we have made use of Eq. ©. Alternatively, we may first perform the 
matching in a non explicitly gauge-invariant fashion, as it is customary in perturbative cal¬ 
culations of Green functions, and then impose gauge invariance at the level of the pNRQCD 
Lagrangian through held redefinitions of the type ©-©• 

Equation (EBP may be generalized to all orders as 
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cr(rXrXB) 


r 2 K 


( 0 )/ 


(70) 


where Id) 0 ' 1 ' stands for dV^/dr. Also Eq. (I69D is valid to all orders. The proof proceeds as 
follows. 

(1) The matching can be performed order by order in 1/m. NRQCD amplitudes that 
may contribute to the matching involve either insertions of two of the operators D 2 /2 m, 
Cp <x-gB/2m and c e ™ cr-eegB cm /2m, or one insertion of the operator — c s a- [— zDx, (/E]/8m 2 
on either the quark or the antiquark line. Couplings of the photon to massless quark loops 
or charm-quark loops may be neglected by the same arguments given in the previous section 

cr-B 

for the matching of V s m . 

(2) First, we consider amplitudes with an insertion of the operator cf 11 cr ■ eeQB em /2m 
and one of either D 2 /2m or cf cr ■ gB/2m. In the Erst case, due to the ultrasoft nature of 
the external photon, we may neglect the action of V on B cm . Hence the magnetic dipole 
operator behaves like the identity operator, and, following an argument similar to that one 
developed in the previous section, we can show that there is no contribution to the matching. 
In the second case, the QCD part of the amplitude factorizes in a term proportional to the 
expectation value of the chromomagnetic operator, which vanishes for parity. 

(3) In all the remaining terms, the electromagnetic coupling is embedded in a covariant 
derivative. Since by projecting the operator 


J d 3 x (-iV - ee Q A em ) ^ + [$ -> ia 2 X *, A em -> -A em ], (71) 

onto the state SHD and taking into account gauge invariance (see Eq. P51Q we obtain 5 
[ d 3 Rd 3 r Tr {S 1 (-iV + r x ee Q B cm + 0(r 2 )) S} , 


5 From 


exp J A em • [~iV xl - eegAfau) - iV X2 + ee Q A(x 2 )\ exp 

and multipole expanding one ends up with — zV + r x eegB em + 0(r 2 ). 
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we conclude that the matching coefficients of the operators 

1 


and 


4 m 2 r 3 
1 


4 m 2 r 3 


{Stcr ■ [r x (r x ee^B cm )]} S 


{St,<r}-[rx(-iV)]S 


( 72 ) 


(73) 


are equal. Indeed, the first is obtained from the second by replacing —iV by r x eegB cm . 
The operator -——- {S^er} • [r x (—iV)] S is protected by Poincare invariance (Grornes 
relation) 


M 
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281 ]. Its matching coefficient is equal to r 2 V^' /2 to all orders in perturba¬ 


tion theory and non-perturbatively. This proves Eq. (|751)- It also proves Eq. flSED, because 

a ■ B 

the analysis does not reveal any contribution to V s m ~ . 

Equations m and (IHHll are both valid to all orders in perturbation theory and non- 
perturbatively. The first equation confirms earlier findings in phenomenological models 
(see, for instance, lizj]). The second one states that to all orders in the strong-coupling 
constant and non-perturbatively the existence of a magnetic coupling of the type induced 
by a scalar interaction is excluded. Phenomenological models often assume that the rela¬ 
tivistic Hamiltonian contains a scalar interaction y°( 1 )H scalar (r) 7 °( 2 b The non-relativistic 
reduction of this term generates, among others, a magnetic spin-flipping term of the type 
^scalar jgt^ . eec? B cm } S. Our analysis shows that such a term is excluded from pN- 


2 m 2 


RQCD. We conclude that a scalar interaction would induce a Ml coupling that in QCD 
cannot be generated, even dynamically, for heavy-quark bound state systems. 6 
The correction induced by the operator 


1 rV^'ir) 


Am 2 2 

to the Ml transition amplitude is 


{S t , cr ■ [r x (r x ee Q B cm )]} S 


1 


A [n 3 Si(0, A) —> n /1 S' 0 (-k) 7 (k, cr)] = -—^—(n'S\rV^ 0) '\nS ). 


6 771 


(74) 


(75) 


Finally, we note that amplitudes with one insertion of the operator —c| m cr ■ [—iV 

cr • V x E cr-VrXr-VS 

x, eegE em ]/8m 2 only contribute to the matching of m2 and V s rn . Since the op¬ 
erator factorizes, there are no soft contributions to the matching coefficients and Eqs. m 
and m turn out to be valid to all orders in a s . 


The situation is here different from the case of the spin-dependent potentials. There, a spin-orbit potential 

The particular nature 


of the type induced by a scalar interaction may be dynamically generated s u 
of the scalar interaction contribution to Ml transitions has also been discussed in |31|- 
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D. Comment on the matching in the strong-coupling regime 


In the weak-coupling regime, at relative order v 2 , the only 1/m 3 operator relevant for Ml 
transitions is' 

jAj {su « 0 B»}v?s. 

Note that corrections to the matching coefficient are suppressed by powers of a s (l/r) <C 1. 
Comparing our expression of the pNRQCD Lagrangian (FTTH) with the phenomenological 


Hamiltonian used in 


12j . we observe that, up to the scalar interaction term, the two expres¬ 


sions are equal. The absence of a scalar interaction in pNRQCD has been discussed above. 
Here, we remark that our expression is valid in the weak-coupling regime (mv 2 > Aqcd) 
only, while phenomenological Hamiltonians are supposed to be applicable to both weakly 
and strongly coupled quarkonia. We may ask how the pNRQCD Lagrangian would change 



liis has been discussed in the absence of an 
35(. Here, we focus on the magnetic dipole 


in the strong-coupling regime (mv ~ A, 
electromagnetic interaction in [30), |32j, |33j, |3; 
couplings. We have shown that the 1/m and 1/m 2 matching is valid beyond perturbation 
theory. However, this is unlikely to happen at order 1/m 3 . Since a s (l jr) ~ 1 is no longer 
a suppression factor, we expect more NRQCD amplitudes to contribute to the matching. 
Among them, we may have amplitudes made of two insertions of the operator cp e r ■ gB/2m 
and one of c'p 11 cr ■ eeQB em /2 m, or one of — c s u • [—iDx,gE]/8m 2 and one of D 2 /2m, or 
one of —c s cr ■ [—iDx, gE]/8m 2 and one of Cp cr ■ gB/2m and so on. These amplitudes will 
be encoded in the matching coefficients of pNRQCD in the form of static Wilson loop am¬ 
plitudes with held strength insertions of the same kind as those that appear in the QCD 
potential at order 1/m 2 Also, they may induce new operators in pNRQCD. A non- 

perturbative derivation of the pNRQCD Lagrangian coupled to the electromagnetic held at 
order 1/m 3 has not been worked out. In a purely analytical approach, such a computation 
will likely have a limited phenomenological impact, due to the many non-perturbative pa¬ 
rameters (Wilson-loop amplitudes) needed. However, if supplemented by lattice simulations, 
it will pave the way for a rigorous QCD study of relativistic corrections to Ml transitions 
in excited heavy-quarkonium states. 


7 The matching coefficients of the operators l/(4m 3 ) {S'/ u l eeQB cmj }V).V/S and l/(m 3 r 2 ){S'l', cr • 
eeQB cm }S are suppressed by powers of a B . Therefore, these operators contribute at relative order a s v 2 
or smaller. 
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In summary, phenomenological models used so far to describe magnetic dipole tran¬ 
sitions in quarkonium, once cleaned of the scalar interaction, appear to be valid only for 
weakly-coupled resonances. For strongly-coupled resonances, at order 1/m 3 , more terms and 
matching coefficients with, in principle, large non-perturbative corrections are expected. 


IV. WAVE-FUNCTION CORRECTIONS TO MAGNETIC DIPOLE TRANSI¬ 
TIONS 

Corrections to the wave function that give contributions of relative order v 2 to the tran¬ 
sition amplitude are of two categories: (A) higher-order potential corrections, which may be 
further distinguished in (A.l) zero-recoil corrections and (A.2) recoil effects of the final-state 
quarkonium, and (B) higher Fock-state corrections. 

A. Corrections to the wave function from higher-order potentials 

1. Zero-recoil effects 

We first consider corrections coming from higher-order potentials that do not depend on 
the centre-of-mass momentum of the recoiling quarkonium. Since 8Vg is, at least, of order 
u 4 , we only need to take into account the correction induced by Eq. E2) to the leading 
amplitude computed in Sec. Ill El The correction is proportional to 

(1 - 6„' n ) (n'5|Tr {{cr, 8V s e n 3 Sl (A) • cr} - 8V S {cr , e n 3 Sl (A) ■ er} jins') . 

It vanishes for allowed magnetic transitions. This follows from the fact that cr ■ eeQB cm /m 
is independent of r and the first-order correction is orthogonal to the zeroth-order wave 
function (see Eq. (EZD). The two terms in the trace come from the correction to the in¬ 
coming and outcoming heavy quarkonium respectively. We distinguish different cases, (i) 
If 8Vg is spin independent, then the trace vanishes, (ii) If 8V$ is a spin-orbit potential, 
then (n'S\8Vs\nS) ~ (n'S'|r x plnS 1 ) = 0 on S' waves, (in) If 8Vs is a tensor potential, 
then (n'S\8Vs\nS) ~ (n / S'|3r*r 7 — r 2 8 l f nS) = 0 on S' waves, (iv) The only non-vanishing 
contribution comes from the spin-spin potential: 

U ss irl 

Tr{S t hU 5 S| = -W TrjSVSo-n . (76) 

i. j 4 m 2 1. J 
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It induces the following correction to the transition amplitude: 


A [n 3 Si(0,A) n' 1 S„(-k) 7 (k,<T)] = c* F "(l - fc„) • ( 77 ) 

m En E n , 

The correction is only relevant for hindered Ml transitions and, in this case, is of order v 2 
(V ss /(Ei 0) - E^) ~ mv 4 /mv 2 ~ v 2 ). 

2. Final-state recoil effects 

The final-state quarkonium is not at rest. It moves with a velocity — k with respect to the 
centre-of-mass frame. In Q, it has been pointed out that due to this motion, higher-order 
potentials that depend on the centre-of-mass momentum may modify the wave function of 
the recoiling quarkonium such that the El operator may induce an effective Ml transition. 
The leading potential relevant to our case is: 

1 

Tt {S f <5V s S} = Tt {{S 7 , <r} • [f x (-tV)] S} . (78) 

We have discussed the spin-orbit potential in Sec. IIII Cl (see Eq. (El), where we noticed 
that its value is protected by Poincare invariance. Inserting Eq. © into Eqs. © and © 
we obtain 

|JZ(P.A»« = -/*fl/<f 3 re iP ' R 'Iv{AA{St,<r} * (W*H(A))|0)} , (79) 

where we have used 

(0, 1 ( 0 ) <ff'(A')|rVf'lff(A)) = * (ff'(A')|p|i/(A)>, 

Eh E h , 

which follows from [p, /i®] = —irV^'^r), and 

E 4>g!(V)( r ) ! (-ff'(A')lplff(A)) = V r *g} A) (r), 

H'jtH, X 1 

which follows from completeness and the definite parity of the functions . Two different 
derivations of Eq. ©, one that uses Lorentz-boost transformations and another one based 
on relativistically covariant formulations, can be found in Appendix IT1 

Equation © states that, due to the recoil, the final state develops a nonzero P- wave, 
spin-flipped component suppressed by a factor vk^/m. As a consequence, in a n 3 Si —> 


25 








?i /1 S'o7 transition, the P -wave spin-triplet final state component can be reached from the 
initial 3 S\ state through an El transition, mediated by the operator 

= j d 3 r Tr {S f r ■ ee 0 E em S} . (80) 

Since the El operator is enhanced by 1/v relative to the leading Ml operator (IH-Hh the 
recoil correction is of order k 1 /m with respect to the leading term. At relative order v 2 , this 
correction is negligible for Ml allowed transitions (fc 7 -C mr 2 ), but should be considered for 
Ml hindered transitions, where ~ mv 2 . 8 

The correction to the transition amplitude is given by 

- W(n ,1 5'o(“k) 7 (k,a)| J d 3 RC^ |n 3 ^(0,A)), (81) 

where |?r' 1 S'o(—k))^-* can be inferred from Eq. CD. After a straightforward calculation, we 
obtain 

A [n 3 Si(0, A) -> n /1 S' 0 (—k) 7 (k, tr)] = (^ n ' n + l - (n'S\r • p|/rS')^ . (82) 

B. Color-octet effects 

In Sec. ITTD21 we pointed out that a heavy quarkonium state also contains higher Fock 
components, in particular, components made of a quark-antiquark pair in an octet con¬ 
figuration. Color-octet effects are regarded as one of the most distinctive benchmarks of 
NRQCD, and have been found to play a crucial role in several phenomenological applica¬ 
tions, e.g. heavy quarkonium decays and productions P- Indeed, color-octet effects are not 
included in any potential-model formulation and have not been considered so far in radiative 
transitions. A color-singlet quarkonium may develop a color-octet component by emitting 
and reabsorbing an ultrasoft gluon. A Ml transition may occur either in the color-singlet or 
in the color-octet component. If Aqcd ~ mv 2 , the process involving a singlet-octet-singlet 
transition is suppressed only by a factor v 2 with respect to the leading one, and, therefore, 
relevant to our analysis. 

8 In the gauge non-invariant formulation of Appendix El the leading El operator is given by Eq. (IBdll . At 
relative order v 2 , it contributes both to allowed and hindered transitions. 
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FIG. 5: Octet contributions to Ml transitions. The single and double lines represent the singlet 
and octet fields, respectively. The circled cross stands for the vertex induced by the interaction 

(Hj . 

In 0,0|, the effect of octet components to the spectrum has been throughly investigated. 
The leading effect is given by 9 


6E„ W = j J dt (0|jE°(R,0)^(0, t)^gE\R, «)|0) r|//(A)>, (83) 

where 0(0. t)^ is a Wilson line in the adjoint representation connecting the point (R, 0) to 
(R, t). Note the appearance of the non-local condensate (0|r/E“(R, 0)0(0, £)^</E fe (R, t) 1 0) , 
typical of the situation Aqcd ~ mv 2 . From Eq. (USD and (jHSD we may calculate the state 
normalization factor 5Zh'- 


_ A ( C 1 (®) Z, (0) \ + 


firr i roc 

sz mx) = —= 7 / dt t (0|jE“(R, 0)0(0, ftfc/E' , (R, f)|0) 

oE k h ’ 0 j 0 


x(//(A)|re- i(£ »’-'“o’ |, r |//(A)>. 


(84) 


The leading color-octet contribution is induced by the chromo-El operator (2BD- At 
second order in r, there are four diagrams contributing to the transition amplitude, as 


Note that in a different normalization for the state |0) is used: |0) = |0)|^ ggj/yTVc. 
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shown in Fig. 0 The contribution of Fig. EJa) corresponds to 


- ( 0 ) < n n S 0 (P')'Y(k,a)\J d^Sil^O, A))j 2) 

- | | 2) (?r ,1 S' 0 (P / ) 7 (k, a) | J d 3 RC ( S 1 \n 3 S l (0, X))^ 

= (2tt) 3 5 3 (P / + k) M (0) [n 3 5!(0, A) n' 1 5 0 (-k) 7 (k, a)] 

A /*oo 

xcT^f J dtt(0\gE a (R, 0)0(0, tfih,E fe (R, t)|0> 
x(7i'S\re~ l< ' E " r \nS) . (85) 


This diagram only contributes to Ml allowed transitions. The contribution of Fig. Mb) 
corresponds to 


-W(n ,1 5 0 (P')7(k,ff)| j d 3 i?££ ) i OCtet |n 3 ,Si(0, A))« 

= (27r) 3 5 3 (P ; + k) M (0) [n 3 S' 1 (0, A) -> n /1 S'o(-k) 7 (k, a)] 


XCri 11 - 


i r°° e H!$>-Ep})t _ x 


dt 0) (0|t/E a (R, 0)0(0, OE*(R, t)|0) 

T/n — EX 


xin'Slre.-^-bo^rlnS) , 

where C^ctet = J d 3 r Tr j^- {0\ a ■ ee Q B em } o|. 

The sum of the contributions of Figs.0D(c) arid 0A1) gives: 

—^(?r' 1 S' 0 (P / ) 7 (k, a) | J ^RC^S^O, A))f 

-5 L 2) (n' 1 ^ 0 (P')7(k,a)| J d 3 R C^S^O, A))< 0) 

= (27r) 3 d 3 (P' + k) M (0) [n 3 <Si(0, A) ?z' 1 S' 0 (-k) 7 (k, a)] 


XlCr 


1 -S n 


/ dt 0) ((%E%R, 0)0(0, t)$g E b (R, t) |0) 

' 0 h/n / 


x(n / S'|r (e 




*o >- — e ) rlnS 1 ). 


( 86 ) 


(87) 


These diagrams only contribute to Ml hindered transitions. 

Remarkably the sum of all octet contributions at relative order v 2 , i.e. Eqs. (ESD-iZD. 
vanishes. This relies, again, on the fact that the leading Ml operator behaves like the identity 

cr-B 

operator in coordinate space. Our cancellation proof, indeed, very much resembles the V s m 
non-renormalization proof given in Sec. 1111 Bl The cancellation of all octet contributions at 
relative order v 2 is nothing more than a manifestation of the proper normalization of the 








quarkonium state when higher-order Fock-space components are taken into account. 10 

Finally, we note that, since the leading operator responsible for El transitions, given 
in Eq. m, is not a unit operator in coordinate space, the above cancellation mechanism 
does not apply there. In the El case, when Aq CD ~ mv 2 , we will, in general, expect octet 
corrections of the same size as the leading relativistic ones. These have not been calculated 
so far, but may be of relevant phenomenological impact for electric-dipole transitions in 
weakly-coupled quarkonia. 

V. FINAL FORMULAE INCORPORATING ALL 0(v 2 ) CORRECTIONS 

In this section, we sum all previously calculated contributions for transition amplitudes 
and give expressions for the widths. 


A. 0(y 2 ) 5-wave transition amplitude 

Summing Eqs. (1351) . (1571) . (1551) . (1551) . (15T1) . (155j) . (1551) . (1751) . (1771) . (1551) . we obtain the Ml 

transition amplitude with 0{y 2 ) corrections included: 


A [u 3 5i(0, A) ->• n /1 S’ 0 (-k) 7 (k, a)] 

— finri (1 + ^Q 1 ) 

l + n e Q (n'S\V ss (r)\nS) 


+(i — 8 nn >)- 


m 


pi(0) p(0) 

- hj n , 


in'S | 


1 4 - 

1 + 1.2 '2 


24 


kT ~ 1 3 +kQ 


„2 if®n 

p + ^-rvp' 


| nS ). (88) 


We have used c| m = 1 + 2 Kg 11 and 


• k / /, | v / n l oP , rV * 1 I \ , tnt ^ 
i—(n r ■ p n = (n —2—- -|- 1 l " x 1 ' 

m \ m z m 


do)/' 


| n) + 0(v A ), 


(89) 


which holds for (1 — 8 nn >) k ~ (1 — 5 nn i) ( En ^ — E®)) ~ (1 — 5 nn i) mv 2 and 5 nn > k ~ 5 nn ’ mv 4 , 
which we neglect. Terms of the type S nn / k, which are beyond our accuracy, have been 
neglected also in Eq. 


10 On this see also Ref. 371. 
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B. S -wave transition widths 


Inserting (18811 into m, we obtain the total transition width for magnetic dipole transi¬ 
tions of S-wave quarkonium np to order k 3 v 2 /m 2 . We use the explicit one-loop value of Kq 11 
given in Eq. ©. In accordance to the power counting, we neglect order a 2 (m) and order 
a s (m)v 2 corrections. The final result reads 


r — 1 2 hj- 

t n 3 Si—7 2 

3 v 77l z 

4 2 k 3 


n 3 Si^n' 1 S 0 'y „ ^ „ 


' Til 


1 + C' 


(n'S I 


«s (m) 


7T 

k 2 r 2 


3 


5 p 2 

6 m 2 


| nS) 


1 (n'S|E ss (r)|mS) 
m 2 

rj n , 


(90) 


(91) 


For completeness, we also give the t^Sq —*► n' 3 S\ 7 transition width, which is relevant only 
for hindered transitions: 


k 3 

r r i 1 5o^n' 3 5i7 n ^ n '4oieQ— ( n ^I 


fc 2 r 2 5 p 2 
7 \nS) 


24 6 m 2 


1 (nW|i/ ss (r)|mS) 

771 2 /A 0 ) _ 

(92) 

Equations very much resemble those derived in Q| and subsequently used 

in most non-relativistic potential model calculations of the magnetic dipole transitions in 
quarkonium (see, for instance, the review in 21]). There are, however, some differences that 
we have already mentioned, but we would like to stress again. 

(1) Equations m-m have a limited range of validity that the EFT framework clarifies. 
They are valid only in the weak-coupling limit, i.e. for quarkonia that fulfill the criterion 
mv 2 > Aqcd- The lowest bottomonium states and the charmonium ground state may belong 
to quarkonia of this kind. As discussed in Sec. IIII Dl in the strong-coupling regime, i.e. for 
higher-quarkonium excitations, at order k 3 v 2 /m 2 more terms will, in principle, arise. 

(2) Equations (190(1 - (192(1 do not contain contributions from a scalar interaction (propor¬ 
tional to — (n / S'|E scalar /m|nS')). This has been often used in potential models, but the 
analysis of Sec. IIII Cl has excluded such a contribution in pNRQCD. 
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(3) The analysis in Sec. lIII Bl has also excluded (to all orders) contributions to the quarko- 
nium magnetic moment coming from the soft scale. This allows us to substitute Kq 11 with the 
value inherited from NRQCD, which at one loop is k < q 1 = Cpa s /The renormalization 
scale of a s is m. 

(4) In Sec. ITYm it has been shown that color-octet contributions, not accessible to 
potential model analyses, cancel at order fc 3 u 2 /m 2 . This leads to the conclusion that in 
the weak-coupling regime, at order k^v 2 /m 2 , Ml transitions are completely accessible to 
perturbation theory. In particular, once the spin-spin potential is written at leading order 
in perturbation theory (which is sufficient here), 

R ss (r) = ^7rC F a s 5 3 (r), (93) 

and Eqs. ed-ed are calculated for Coulomb wave functions, the transition rates will only 
depend on the strong-coupling constant. 


C. P -wave transition widths 


In this section, we consider only allowed Ml transitions between P -wave states, since 
hindered P -wave transitions are unlikely to accommodate within a weakly-coupled picture. 
The calculation proceeds very much like the analogous one for S'-wave states, so we will 
not present details here. Octet contributions again cancel by the same argument as given 
for S'-wave transitions. At order fc 3 v 2 /m 2 , only two operators contribute to Ml allowed 

1 rVjp\r) 


transitions: 


{S t ,cr.ee 0 B cm } V^S and 


{s', 


er • r x 


4m 3 1 ’ ) r 4m 2 2 

Summing their contributions, at order fc 3 r; 2 /m 2 , the final results read 


r„3 P 


_ 4 2 kj- 

'W p l7 — g ae Q m 2 


1 +C F 


«sH j /.mP* 


7T 


-dj ( n P\d-\nP) 
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4 k 3 

ritlp^naPj-, = (2J+l)-£*eQ—^ 


1 + C F ^P - dj (nP\^\r,P) 
5T m 


(txee 0 B'")]}S. 


(94) 

(95) 


where do = 1, d\ — 2 and d 2 = 8/5. We have made use of the virial theorem. Corrections 


induced by the operator 


1 rVt 


d°)'( r ) 

{S f , cr • [f x (r x ee Q B cm )]} S vanish for J = 0 


4m 2 


states. 

Combining Eq. with Eqs. (ED and ED, we obtain that, at leading order, the following 
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relations hold: 


3 r ri ip 1 ^ n 3 Po7 r n 3 5l _ n i 5o7 

^n^P2~>n 1 Pi 7 r n 3 5 l ^ n i 5o7 


= 10 , 


— m 1 *S'o 7 ^n^Pi—^n^Pi 7 

r n3P2 

—miPi 7 r„3 S i 

—>n 1 So 7 


= 5, 


(96) 


P“ p“ 

which follow from (n5|—-|n5) = (nP —-| nP) = 


E< 


(o) 


11 


771 - 


771- 


771 


VI. APPLICATIONS 

We have remarked that Eqs. m-m, © and (USD are valid only for weakly-coupled 
quarkonia. It is generally believed that the lowest-lying bb states, T(IS') and r)b, are in the 
weak-coupling regime. The situation for x&(lP), hb(lP), T(2 S) and 7]b(2S) is more con¬ 
troversial, as it is for the lowest-lying cc states. We will assume that also these states are 
weakly coupled and see whether the comparison between our predictions and the experi¬ 
mental data supports this assumption or not. As for the n — 2 charmonium states, it is 
undoubtedly inappropriate to consider them as weakly-coupled systems. A further compli¬ 
cation of the t(j(2S) and r] c (2S) states is that they he too close to the open charm threshold, 
so that threshold effects should be included in a proper EFT treatment. We will not con¬ 
sider them in our analysis. In the following, we shall apply Eqs. (FIHll - dU^ll to J/%b —■> 7] c 7 , 
T(IS') —■> r/fc 7 , T(25) —■> 76 ( 25 ) 7 , T(25) —■> r/b 7 and rjb(2S) —> 7(15)7, and Eqs. (liTEl) and 
© to h b (lP) -> X6o,i( 1-P)7 and Xb 2 ^P) -> hb(lP) 7 . 

A. J/V’-^7c 7 


In potential models, the transition J/-0 —> rj c 7 has been often considered problematic to 
accommodate because its leading-order width is about 2.83 keV (for m c = Mj/^/2 = 1548 
MeV), far away from the experimental value of (1.18 ± 0.36) keV 0 . 

Since we assume that the charmonium ground state is a weakly-coupled quarkonium, 
Eq. m provides the transition width up to order k^v^/ 771 2 . We may conveniently rewrite 


11 So far, we have labeled P -wave states with their principal quantum number n. In the next section, we 
will follow the usual convention for which a y(!P) state is a n = 2, L = 1 state. 


32 







it as 


r 


_ 16 2 k* 

J/ip^v c 7 — g ae ° M 2 




16 2 

= —o / p z -- 

3 C M 2 


J/i/j L 


' g g (M w /2) 2 (lS\3V^ ) ~rV^ , \lS) 
F vr 3 Mj/^ 

1 + C F a ‘ {Mj, * /2) - \(C Fa ,( P]/il )f 

7r o 


( 97 ) 


where in the first line we have reexpressed the charm mass in terms of the J/ip mass, 


P z 


M m = 2 m c + {1S\^~ + V<°\r)\lS), 

TYl r 


and made use of the virial theorem to get rid of the kinetic energy. We have made explicit 
in Eq. m that the normalization scale for the ck s inherited from is the charm mass 
(a s (Mj/y/2) as 0.35), and for the a s , which comes from the Coulomb potential, is the typical 
momentum transfer pj/^ ~ mCFOt s {pj/ip )/2 ~ 0.8 GeV. Numerically we obtain: 


r. w ^„ o7 = (1.5 ± 1.0) keV. 


(98) 


The uncertainty has been estimated by assuming the next corrections to be suppressed by 
a factor oi^(pj/^) with respect to the transition width in the non-relativistic limit. 

Some comments are in order. First, we note that the uncertainty in (jHHl) is large. In 
our view, it fully accounts for the large uncertainty coming from higher-order relativistic 
corrections, which may be large if we consider that those of order k^v 2 /m 2 have reduced 
the leading-order result by about 50%, and for the uncertainties in the normalization scales 
of the strong-coupling constant. Both uncertainties may be only reduced by higher-order 
calculations. 

Despite the uncertainties, the value given in Eq. (10811 is perfectly consistent with the 
experimental one. This means that assuming the ground-state charmonium to be a weakly- 
coupled system leads to relativistic corrections to the transition width of the right sign and 
size. This is not trivial. If we look at the expression after the first equality in Eq. (EH) , we 
may notice that 3V^°' ) — rV^' is negative in the case of a Coulomb potential (i.e. it lowers 
the transition width), but positive in the case of a confining linear potential (i.e. it increases 
the transition width). This may explain some of the difficulties met by potential models in 
reproducing Tj/^ Vc , r In any rate, it should be remembered that Eq. (l97l) is not the correct 
expression to be used in the strong-coupling regime. 
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FIG. 6: r T (is)^r /( , 7 and r T ( 2 S)^%( 2 S )7 as a function of the photon energy. 

B. T(lS)-7fe 7 , T(2S)^ Vb (2S)^ 


Allowed Ml transitions in the bottomonium system that may be treated by the weak- 
coupling formula m are —>■ Vb 7 and, perhaps, T(2S) —> tjb(2S) 7 . We have 


16 


P 

2 K i 


r T(lS)-»Tft 7 - M'- 
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T(2S)^%(25)7 


= — ae 


T(1S) 
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1 + Cp « i pWH)_2 , ' 
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(99) 


1 + Cj 


a s (Mr(is)/2) f Cpcxs(pr(is)) 


7T 


5 / C , f« s (pt(25)) \ 2 

3 V 4 J 


( 100 ) 


where we have expressed the b mass in terms of the T(1S) mass. We have made explicit that 
the renormalization scale for the a s , inherited from K b m , is the bottom mass (« s (Mx(is)/ 2 ) ~ 
0 . 22 ), while for the a s , which comes from the Coulomb potential in the T(IS') system, is 
the typical momentum transfer pr(is) ~ mCFOL s (p-t{is))/2 ~ 1.2 GeV, and for the a s , which 
comes from the Coulomb potential in the T(2S) system, is the typical momentum transfer 
p x(2S ) ps mC F a s (p T( 2S))/4 ps 0.9 GeV. 

Since the 7] b has not been discovered yet, in Fig. El we show rx(is )_> %7 and Tx^sp^psp 
as a function of fc 7 . The bands stand for the uncertainties calculated as the product of the 
transition widths in the non-relativistic 
we use the value of the r) b mass given in 


irnit by p (prps)) an d pUPps)) respectively. If 


3 


i.e. 


k~ = 39 ± 13 MeV, we have 


rT(isp %7 — (3-6 ± 2.9) eV, 


( 101 ) 
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FIG. 7: r t( 2 S)— 7 and r, % ( 2 s)^x(i 5)7 as a function of the photon energy, 
which corresponds to a branching fraction of ( 6.8 ± 5.5) x 10~ 5 . 


C. T(2S)-, %7 , % (2S)->T(1S) 7 


For hindered Ml transitions, Eqs. m and m only provide the leading-order expres¬ 
sions. We consider, here, Y (2 S) —> rjb'y and rjb(2S) —> Y( 1 S ')7 transition widths that we 
write as 


rx(25)^ %7 

r % (2S)^T(lS)7 





M 2 

lvl r(is) 


16a?e 2 


k 3 

^T(lS) 


3lV2 2 1024V2 k 2 

81 1 Fasj + 729 (M X(15) C F a s )2_ 

41vg 2 1024^2 k* 

81 1 ™ + 729 (. M r{ is)C F a s ) 2 


( 102 ) 

(103) 


Since terms arising from the IS 1 and 2S system mix, it is difficult to assign a natural normal¬ 
ization scale to a s appearing in Eqs. (EH) and 'EE' without doing a higher-order calculation. 
In Fig. d we show a plot of rx( 2 S )-> % 7 and r % ( 2 5 )^T(is )7 as a function of fc 7 . The scale of 
ct s appearing in Eqs. EH and EH has been arbitrarily fixed to 1 GeV. The bands stand 
for the uncertainties calculated as the products of the transition widths by a? s (px( 2 S))- 
CLEO III recently set the 90% upper limit for the branching fraction of T(2S) —s► 7^7 
to be 0.5xl0 ” 3 0 The values plotted in Fig. dare about a factor 10 above the limit . 12 


12 


Large contributions stem from the spin-spin potential term. If instead of using E^ 
we use the physical mass difference, the decay width reduces by about a factor one half. 


in this term, 
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FIG. 8 : r /lb( 1 p)^ Xb0 ( 1P ) 7 , r^ ( 1 p)_ X 6 l( 1 p ) 7 and ^ Xb2 {iP)^ hb (iP) ^ as a function of the photon energy. 
The vertical lines correspond to the centre-of-gravity mass of the \bj(lP) states, which is about 
9900 MeV. This is believed to be a rather accurate estimate of the h b (lP) mass. For this value 
of the h b {lP) mass, we obtain ^ hb (iP)^ Xb0 ^P)i = 1 ± 0.2 keV, T hb(1P) ^ Xbl{1P ) ^ = 17 ± 4 eV and 
T Xb 2 (iP)^h b (iP)'y = 90 ± 20 eV. 



Despite the fact that our calculation is just a leading order one and, therefore, potentially 
affected by large uncertainties, it is not obvious that perturbation theory may accommodate 
for such a large discrepancy. In case, this may hint to a strongly-coupled interpretation of 
the bottomonium 2 S states. 13 14 


13 A conclusion of this kind has been reached in from the study of T (2 S) —> A' 7 radiative decays. 

On the other hand, the masses of the n = 2 bottomonium states seem easier to accommodate within a 
weakly-coupled picture Q| ■ 

14 It has been also argued that new physics may broaden the r]b resonance, which thereby may have escaped 
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D. h b (lP) ->• Xbo,i( lp ) h Xb2(lP) -> /ift(l-P) 7 


P-wave Ml transitions that may be possibly described by pNRQCD in the weak-coupling 
regime are Ml allowed transitions between n — 2 bottomoninm states. Proceeding like in 
Sec. IV1BI we obtain 


M1P)^X6o(1P)7 
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and 


16 2 k i 

V Xb2iXP)^h b {iP) 1 - ~^-ote b 
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° JW T(15) 
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(106) 


Since the h b (lP) has not been discovered yet, in Fig. |S] we show r/ l6 ( 1 p)_ >X60 ( 1 p) 7 , 
r^ ( ip)^ XM( ip )7 and r XM (ip)_ h 6 (ip )7 as a function of fc 7 . We assume Px( Xi ,(iP)) ~ Pt(t(2S)) ~ 
0.9 GeV. The bands stand for the uncertainties calculated as in Sec. IVI Bl Specific predic¬ 
tions of pNRQCD in the weak-coupling regime are also Eqs. © for n — 2. 


VII. CONCLUSION AND OUTLOOK 

The paper constitutes a thorough study of magnetic dipole transitions in the framework 
of non-relativistic EFTs of QCD and, in particular, of pNRQCD. The matching of the 
magnetic dipole operators at order 1/m and 1/m 2 of pNRQCD has been performed to all 
orders in a s . The matching at order 1/m 3 has been carried out at leading order in the 
weak-coupling regime. Relativistic corrections to the transition widths have been included 

detection at CLEO III 0|. 
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in a systematic fashion. Having achieved this, we could answer the questions raised in the 
introduction, (i) The contribution to the quarkonium anomalous magnetic moment coming 
from the soft scale vanishes to all orders, (ii) There are no contributions to the magnetic 
dipole operators of the type induced by a scalar potential. (Hi) In the weak-coupling regime, 


non-perturbative corrections due to color-octet con 


;ri 


Dutions vanish at relative order v 2 . 


Our final formulae ©-© are the same as in |l5|, once cleaned of the scalar potential 
and once the one-loop expression of the quarkonium anomalous magnetic moment has been 
used. They are valid only for quarkonia that fulfill the condition mv 2 > Aqcd, he. only 
for the lowest-lying resonances. The application of Eq. m to the transition J /-0 ^ rj c 7 
shows that a weak-coupling treatment of the charmonium ground state is consistent with 
the data. We also provide a prediction for the analogous transition in the bottomonium 
case. Equations © and © are, to our knowledge, new. 

Higher resonances that obey the condition Aq CD ~ mv are described by pNRQCD in the 
strong-coupling regime. In this case, more operators, arising from the 1/m 3 matching, are 
likely needed. 

This work provides a first step towards a complete treatment of quarkonium radiative 
transitions in the framework of non-relativistic EFTs of QCD. Some of the next steps are 
obvious and we shall conclude by commenting on some of them. 

(1) To describe Ml transitions for higher resonances the completion of the non- 
perturbative matching of the relevant pNRQCD operators at order 1/m 3 will be needed. 
The matching coefficients will be Wilson-loop amplitudes similar to those that describe the 
non-perturbative potential at order 1/m 2 |2£j. This calculation, combined with a lattice 
simulation of the Wilson-loop amplitudes, may provide a rigorous QCD derivation for all 
quarkonium Ml transitions below threshold. 

(2) Ml hindered transitions of the type T(3S') —> ^7 have been also studied at CLEO 

HI |4£j. They involve emitted photons whose momentum is comparable with the typical 
momentum-transfer in the bound state. Then one cannot rely anymore on the multipole 
expansion of the external electromagnetic fields. In this case, however, one may exploit the 
hierarchy Aq CD > pr( 3 S)- T° our knowledge, this situation has not been analyzed so 

far. 

(3) Electric dipole transitions have been mentioned only superficially in the paper. The 
pNRQCD operators, relevant for El transitions beyond leading order, have not been given. 


38 





In the weak-coupling regime, octet contributions may be important and can be worked out 
along the lines discussed here. However, most of El transitions may need to be treated in a 
strong-coupling framework. 
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APPENDIX A: FINAL-STATE RECOIL EFFECTS 


We present here two alternative derivations of the final-state recoil effects calculated in 

Sec.UVX! 

1. Recoil effects from Lorentz boosts 


The effect on the quarkonium state of higher-order potentials that depend on the centre- 
of-mass momentum P may also be calculated by boosting the quarkonium state at rest by 
—P /Mh ~ —P/(2 m) (i.e. minus the recoiling velocity). The importance of boost effects 


12j |. In our 


on the final-state quarkonium was first pointed out by Grotch and Sebastian in 
language, their argument goes as follows. 

The Lorentz-boost generators K of pNRQCD may be read from j3|. The leading spin¬ 
flipping contribution to K is given by 


6K= d 3 R / d 3 r —Tr { [S f , «r x V r ] S} 


It boosts the held St by an amount 
hS f = -i 

Substituting (E2> into 


•hK,S t 
2 m 


£ijk 


8m 2 


{vjsvq 


d?R / d 3 re ipR Tr{* ffW (r)<SS , (r,R)|())} , 


we 


obtain Eq. GSD- 


(Al) 


(A2) 


(A3) 


2. Covariant formulation 


Final-state recoil effects are automatically included in any Lorentz covariant definition of 
the wave function, like that one provided by the Bethe-Salpeter equation |46] . In momentum 
space, the Bethe-Salpeter wave function has the following spin structure: 


, BS If /2 + /> + m 1 Pfi n If/2 — /> — m 




Qh 


(A4) 


2 m 2 2 m 

where, at the order we are interested in, P^ « (2m, P) is the centre-of-mass momentum, 
~ (0, p) is the quark-antiquark relative momentum, « (l,P/(2m)), G^So = 7 5 an< i 
Gn 3Si(A) = fasA A), with u ■ e n 3 Sl {\) = 0. 
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Expanding in P and p and keeping the upper-right 2 <8) 2 block we obtain 



(AS) 



oc <r • e — 


- • (e„3 Sl (A) x p) + ... . 


(A6) 


The Erst terms in the equations give the spin structures of Eqs. (EDI) and m , the second ones 
provide (0|S(r, R)|n 1 5'o(P)) (1 ' ) and (0|5(r, R)|n 3 S'i(P. A)) (1 ) respectively, where \H(P. X))^ 
has been given in Eq. Cl. 

APPENDIX B: GAUGE INVARIANCE 

In the main text, we have employed an explicitly gauge-invariant formulation. In the 
literature, however, this has never been the case. As a consequence, partial results may 
differ. In this appendix, in order to make contact with the existing literature, we recalculate 
Ml transitions in a formulation of pNRQCD where U( l) em gauge invariance is not manifest 
at the Lagrangian level. This means that we shall express the pNRQCD Lagrangian in terms 
of the fields S' and O' defined in Eq. (®. Of course, the final, total results are identical in 
the two formulations. 

If the calculation of Ml transitions in pNRQCD is performed in terms of the held S', 
there are two corresponding changes. 

(1) The first change concerns 1/m 2 operators. As discussed in Sec. IIII 01 these may be 
obtained by projecting CD onto a two-quark state. If the projection is performed on (1551) 
and flU, we obtain the operator 15 



(Bl) 


It induces the following correction to S'-wave transition widths 


A [n s S,(0, A) -» n ,1 S 0 (-k)'y(k, <?)] = -A— {n'S|rUj 0) '|iiS), 


(B2) 


which differs by a factor 1/2 from Eq. (1751) . 


15 The leading operator in the multipole expansion, proportional to <x • r x eeQA em , does not contribute to 
Ml transitions. 
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(2) The second change concerns final-state recoil effects. These have been calculated in 
a gauge invariant formulation in Sec. IIV A 21 In terms of the fields S', El transitions are 
mediated by (to be compared with Eq. (IHU1B 


-2 i / d 3 r Tr S ,f 


, /t ee Q A cm • V,. g , 


m 


(B3) 


The correction to S -wave transition widths induced by dm on a recoiling final state i 


is 


A A) —> n /1 S' 0 (-k) 7 (k, a)] = ~{n S\-^^\nS). 


(B4) 


12]. Note that, at order u 1 2 , Eq. (IB4J) also contributes 


This is exactly the result first derived in 
to Ml allowed transitions, while Eq. (PI) only contributes to Ml hindered transitions. 
Summing Eqs. (B and (IBID we obtain 


(n'S | 


12 m' ' s 

Summing Eqs. m and P|) we obtain 




6m 2 


\nS). 


nS). 

By using Eq. © and S n > n k ~ 5 n > n mv 4 5 6 7 , one can easily see that the two expressions are 
equal at order v 2 . It is straightforward to perform the same check also in the case of P -wave 
transitions. 


(n’S | 


1 T 7-(°) ' , ^ 

--—rVg + -— 
om 4m 


Sn'n + 2 r • P 
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